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Testing new strategies in finite density 

V. Azcoiti a , G. Di Carlo b , A. Galante bc *, V. Laliena a . 
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A new approach for non zero chemical potential simulations is tested in the Gross-Neveu model for infinite 
flavor number, where the critical line is reconstructed in a large [i/T interval. A comparison with results from 
standard imaginary chemical potential approach as well as first results for Nf — 4 QCD are presented. 



Recently we proposed a new method for sim- 
ulate systems at non zero chemical potential pQ. 
The starting point is to note that the temporal 
part of the Dirac operator (the one that con- 
tains the /i dependence) can be written as the 
sum of two contributions: the first, proportional 
to cosh(/z), multiplies an antihermitian operator, 
the second, proportional to sinh(/x), multiplies an 
hermitian operator. Since the spacial part of the 
fermionic operator is antihermitian too, we see 
that for non zero fi the Dirac matrix looses its 
global antihermiticity and this is the origin of the 
complex action problem. 

Clearly the imaginary chemical potential solves 
the problem: sinh(/i) gets an extra factor i and 
converts the corresponding operator to an antier- 
mitian one. Introducing an imaginary chemical 
potential we look for a critical line in the (/3, i/i) 
plane to analytically continue in the (/3, /x) one: 
in doing that we are forced to consider larger val- 
ues of the physical temperature T than the ones 
of physical relevance. 

In our proposal we promote cosh(/i) and 
sinh(/x) to independent variables x and y and ex- 
plore the phase structure as a function of such 
variables to recover the physical critical point 
along the line that corresponds to x = cosh(/i), 
y = sinh(/i). To do that we perform standard 
(i.e. with real and positive action) simulations at 
pure imaginary y and extend the results to real 
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y. This can be achieved in two different ways |T]. 

The first is to fix the inverse gauge coupling 
and (i) find the critical line in the (x, iy) plane, 
(ii) fit this line with an even polynomial of y, (Hi) 
analytically continue the critical line to the (x, y) 
plane, (iv) from the intersection of the analyti- 
cally continued critical line with the x 2 — y 2 = 1 
line (which corresponds to real and positive values 
of fj.) it is possible to read the critical value of /j,. 
Working at fixed lattice temporal extent (N t ) this 
is equivalent to do simulations at fixed T i.e. the 
same physical temperature of the physical critical 
point. 

The second approach is to choose a value for 
[i and fix the x variable to cosh(/x). By standard 
simulations we find the critical line in the ((3, iy) 
plane and proceed in close analogy with previ- 
ous case. Only point (iv) changes: to determine 
the critical value of (3 at given \i we look for the 
intersection of the analytically continued critical 
line with y — sinh(/i). In this case, for a given 
N t , this is equivalent to search for the transition 
point doing simulations at fixed n/T. 

Comparing the standard imaginary chemical 
potential approach [2] and our approach(es) we 
can note that in both case we have to rely on an 
analytical continuation of the critical line (with 
all the uncontrolled systematics that this proce- 
dure can in principle introduce). 

The simulation cost of our method(s) is bigger 
since we need to explore a two dimensional space 
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Figure 1. G-N model: analytical continuation of 
the critical line from imaginary to real y and de- 
termination of the critical \i. 



(the x — iy or (3 — iy one) to extract the critical 
line. Besides that we can expect to be somehow 
closer to the physical critical point (in one case 
we work at same /?, in the other at the same value 
for the adimcnsional lattice fi) and this could im- 
prove to reduce the overall systematic effects due 
to extrapolation (i.e. reach larger fi/T values). 

To check if the above prejudices are correct or 
not we consider an analytically solvable model: 
the Gross-Neveu model (i.e. the four-Fermi model 
in 2 + 1 dimensions) at finite T and /i. 

It can be solved analytically in the large Nf 
limit and its phase structure in the (/x, T) plane 
resembles very much the one expected for QCD, 
with a phase transition line separating a phase 
where chiral symmetry is spontaneously broken 
from a phase where the symmetry is restored 3, . 

After the introduction of an auxiliary scalar 
field a we can define a generalized fermion action 
(we consider the zero mass limit) 
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Figure 2. G-N model: the critical line 



where S S p is the standard spacial part of the 
Kogut-Susskind action and < n,n > denotes the 
set of 8 dual lattice sites n surrounding the site n. 
The above actions agrees with the standard one 
for x = cosh(zx), y — sinh(/i). In the Nf —> oo 
limit the generalized model can be analytically 
solved by the gap equation 



-g 2 (W) 



(i) 



where £ is the expectation value of the scalar 
field. 

We used Q in the infinite spacial volume limit 
to find the critical line in the (x, iy) and (j3, iy) 
planes and then apply our procedures to extract 
the critical line as a function of y. These re- 
sults can be directly compared with the exact 
ones we can extract from the gap equation in or- 
der to check, in a non trivial case, the potential- 
ities of our approach. The final goal is to recon- 
struct the critical line in the (/3, fx) plane or, using 
adimensional variables, in the (/i/£o,T/£o) one 
where Eo is the infinite volume chiral condensate 
at n = 0. 

For Nt = 4 fig. 1 is an example of how 
the method works. The points are the exact 
gap equation solution for the critical line (both 
for real and imaginary y) and a two parame- 
ter quadratic fit is satisfactory (upper continuous 
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Figure 3. G-N model: comparison of our results 
from (x, iy) and (/3, iy) analysis and exact result 
in the /u/Sq — T/Y.Q plane. 



Figure 4. As in fig. 1: determination of the criti- 
cal fj, for Nf — 4 QCD (the symbols are the crit- 
ical points at imaginary y). 



line). The analytical continuation to real y gives 
the decreasing continuous line and its intersection 
with the x 2 — y 2 = 1 line allows to determine the 
critical point. Fig. 2 summarizes the results in 
the (/?, /i) plane compared with the exact ones as 
well as with the imaginary chemical potential re- 
sults based on a quadratic approximation of the 
critical line. 

We get similar results in the (/3, y) plane: the 
results are reported in fig. 3 and compared with 
the (x,y) approach for the N t = 8 case . Results 
from (x, iy) plane are closer to the exact solu- 
tion respect to (/3, iy) ones and, when comparing 
with the standard imaginary chemical potential 
approach (considering in both cases an analytical 
continuation based on a quadratic two parameter 
fit), we see a clear improvement when the critical 
chemical potential is large (see fig. 2). Indeed 
we pick up the correct result up to \ifT ~ 8. 
From considering different Nt values we realize 
that the quality of our results seems to decrease 
increasing Nt and/or increasing the critical value 
of the chemical potential. We considered up to 
Nt = 12 (far beyond any optimistic prevision for 
large scale simulations in finite density QCD in 



the next years) and, up to this value, the quanti- 
tative agreement is good. 

We have the first results for N f = 4 full QCD 
in a 8 3 x 4 lattice at m q = 0.06. The critical 
point is located from peaks in the susceptibility 
of plaquette, Polyakov loop or chiral condensate 
(they all agree). Fig. 4 shows the determination 
of the critical fi based on continuation to real y 
of a quadratic fit to the critical line in the (x, iy) 
plane at j3 = 5.05. A consistent result is obtained 
from an independent analysis carried out from a 
different data set based on numerical simulations 
at fixed n to determine the (/3, iy) critical line. 
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